A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
11.03.5; 11.09.3; 11.05.2; 11.01.1 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



The Influence of Stellar Energetics and Dark Matter on 
the Chemical Evolution of Dwarf Irregulars 

F. Bradamante 1 ' 2 , F. Matteucci 1 ' 2 and A. D'Ercole 3 

1 Dipartimento di Astronomia, Universita di Trieste, Via G.B. Tiepolo 11, 34131 Trieste, Italy 

2 SISSA/ISAS, Via Beirut 2-4, 34014 Trieste, Italy 

3 Osservatorio Astronomico di Bologna, Via Zamboni 33, 44100 Bologna, Italy 
February 1, 2008 



Abstract. A chemical evolution model following the evo- 
lution of the abundances of H, He, C, N, O and Fe for 
dwarf irregular and blue compact galaxies is presented. 
This model takes into account detailed nucleosynthesis 
and computes in detail the rates of supernovae of type 
II and I. The star formation is assumed to have proceeded 
in short but intense bursts. The novelty relative to pre- 
vious models is that the development of a galactic wind 
is studied in detail by taking into account the energy in- 
jected into the interstellar medium (ISM) from both su- 
pernovae and stellar winds from massive stars as well as 
the presence of dark matter halos. Both metal enriched 
and normal winds have been considered. 
Our main conclusions are: i) a substantial amount of dark 
matter (from 1 to 50 times larger than the luminous mat- 
ter) is required in order to avoid the complete destruction 
of such galaxies during strong starbursts, and ii) the en- 
ergy injected by stellar winds and type la supernovae into 
the ISM is negligible relative to the total thermal energy, 
and in particular to the type II supernovae, which in fact, 
dominate the energetics during starbursts. 



1. Introduction 

Dwarf irregulars galaxies (DIG) and blue compact galax- 
ies (BCG) are very interesting objects to study galaxy 
evolution, for they are easier to model because of their 
generally small sizes and simple structures. Their stellar 
populations appear to be mostly young, their metallicity 
is low and their gas content is large. All these features in- 
dicate that these galaxies are poorly evolved objects, and 
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may have undergone discontinuous or gasping star forma- 
tion activity. In particular, BCG are the most luminous 
among dwarf irregular galaxies: their very low metal con- 
tent, large gas fraction and very blue colours, suggest that 
star formation may have proceeded in short and intense 
bursts of activity (Searle et al. 1978). On the other hand, 
DIG seem to show a 'gasping' activity of star formation 
(Tosi et al. 1992). In the last few years a great deal of 
theoretical models for DIG and BCG appeared in the lit- 
erature. By means of analytical models, Matteucci and 
Chiosi (1983) first studied the effects of galactic winds on 
the chemical evolution of these galaxies, and suggested, 
among other hypotheses, that galactic winds of variable 
intensity could explain the observed spread in the physi- 
cal properties of these galaxies. In the following years, by 
means of numerical models, Matteucci and Tosi (1985), 
Pilyugin (1992, 1993) and Marconi et al. (1994), studied 
the effects of galactic winds on the chemical evolution of 
such galaxies and confirmed the previous results. In these 
papers the possibility of metal enriched winds was also 
considered and favoured relative to the normal galactic 
winds. This metal enriched wind hypothesis, in fact, al- 
lowed them to explain the observed He abundances ver- 
sus metal abundances. However, none of these studies had 
taken into account either the energy injected by super- 
novae and stellar winds into the ISM or the presence of 
dark matter halos which are of fundamental importance 
for the development of a galactic wind. The galactic winds 
were just assumed to occur every time the starburst was 
active. 

Dark matter in dwarf irregular galaxies seems to be 
quite important as shown by recent data indicating an in- 
creasing amount of dark matter with decreasing luminos- 
ity (Skillman 1996), and therefore it cannot be ignored in 
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chemical evolution models. Kumai and Tosa (1992) also 
explored the possibility of dark matter in these galaxies 
in a very simple way, and suggested that various amounts 
of dark matter could explain the observed spread in the 
Z versus log(M gas /M tot ) relation. 

The energetics of supernovae and stellar winds and 
their interaction with the ISM, is a difficult problem, still 
not quite understood. The only existing models of chem- 
ical evolution taking into account the stellar energetics 
relative to the occurrence of galactic winds, are those de- 
veloped to study the chemical and photometric evolution 
of elliptical galaxies (Arimoto and Yoshii 1987; Matteucci 
and Tornambe 1987; Angeletti and Giannone 1990; Fer- 
rini and Poggianti 1994; Bressan et al. 1994; Gibson 1994, 
1997). In these models, some specific assumptions rela- 
tive to the amount of energy transferred into the ISM are 
made. Gibson (1994) has shown that the energetics from 
stellar winds is negligible in massive ellipticals whereas 
it could be important in smaller ones. For this reason it 
seems worthwhile to explore the effects of both supernovae 
and stellar winds on the evolution of dwarf irregulars and 
blue compact galaxies. 

The aim of this paper is to take into account both 
the stellar energetics and the presence of dark matter in 
modelling the chemical evolution of dwarf irregular galax- 
ies. The development of a galactic wind will therefore be 
calculated in detail and the existence of metal enriched 
or differential winds will be taken into account. By dif- 
ferential wind we intend a wind carrying out only some 
heavy elements, in particular those produced by type II 
supernovae which are the predominant supernovae during 
a starburst (Marconi et al. 1994) and eject material at 
much higher velocities than normal stellar winds. There- 
fore, under this assumption elements such as nitrogen and 
helium, which are restored by low and intermediate mass 
stars through stellar winds, will not leave the star forming 
region whereas oxygen and the other a-elements, ejected 
during type II supernova explosions are likely to be ejected 
outside the region and perhaps outside the galaxy. Iron 
will only in part leave the star forming region since only 
a fraction of this element is produced by SNe II, whereas 
the bulk of it comes from SNe la. It should be again said 
that the differential wind assumption has proven to be 
the most viable solution to the helium problem in blue 
compact galaxies (see Pilyugin 1993 and Marconi et al. 
1994). 

The paper is organized as follows: in section 2. we will 
present the observational data concerning dwarf irregulars 
and blue compact galaxies, in section 3. we will describe 
the chemical evolution model and the nucleosynthesis pre- 



scriptions, in section 4. we will give our model results and 
finally in section 5. we will draw some conclusions. 

2. Data sample 

The data we refer to are those presented in Marconi et 
al. (1994) (mainly Pagel et al. 1992 and Campbell 1992) 
and in Skillmann (1996) (mainly Thuan 1995 and Garnett 
1990). 

Thuan at al. (1995) presented high quality spectropho- 
tometric observations of 15 supergiant HII regions in 14 
BCG. We compare our model results with their data rela- 
tive to N, O, and Fe. Their data showed that none of the 
ratios N/O, Fc/O, Ne/O, Ar/O, S/O, depends on the oxy- 
gen abundance indicating that abundance ratios in these 
galaxies show the pollution only from massive stars and 
that all the elements, including nitrogen, should have a 
primary origin (see paragraph 3.1). 

Garnett et al. (1995a) presented UV observations of 7 
HII regions in low-luminosity dwarf irregular galaxies and 
the Magellanic Clouds (HST data). This allowed them to 
measure the carbon abundances and discuss the evolu- 
tion of the C/O ratio in these galaxies. They found that 
the C/O ratio increases continuously with increasing O/H 
and suggested that this can be due to the fact that in 
the most metal poor galaxies, massive star nucleosynthesis 
predominates (oxygen production), while in more metal- 
lic objects the delayed contribution of intermediate mass 
stars to carbon is more evident. However, their compari- 
son with chemical evolution models is not really appropri- 
ate since they have adopted, for such comparison, models 
relative to the solar neighbourhood, where the star for- 
mation mechanism and evolution are definitively different 
from those of dwarf irregular galaxies. 

3. Theoretical Prescriptions 

To calculate the chemical evolution of the star forming 
region of DIG and BCG, we used an improved version of 
the model presented in Marconi et al. (1994). The main 
features of our chemical evolution model are the following: 

1. one-zone, with instantaneous and complete mixing of 
gas inside this zone; 

2. no instantaneous recycling approximation, i.e. the stel- 
lar lifetimes are taken into account; 

3. the evolution of several chemical elements (He, C, N, 
O, Fe) due to stellar nucleosynthesis, stellar mass ejec- 
tion, galactic wind powered by supernovae and stellar 
wind and infall of primordial gas, is followed in details. 
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If d is the fractional mass of the element % in the gas, its 
evolution is given by the equations: 



Gi = -mXi{t) + Ri(t) + (Gi)inf - G iw (t) 



(1) 



where Gi(t) = M g (t)Xi(t)/M L (t G ) is the gas mass in the 
form of an element i normalized to a total luminous mass 
M L = 10 9 M Q . 

The quantity Xi(t) — Gi(t)/G(t) represents the abun- 
dance by mass of an element i and by definition the sum- 
mation over all the elements present in the gas mixture 
is equal to unity. The quantity G(t) = M g (t)/ML is the 
total fractional mass of gas. 

The star formation rate we assume during a burst, 
ip(t), is defined as: 



m = TG(t) 



(2) 



where T is the star formation efficiency (expressed in units 
of Gyr~ 1 ) 1 and represents the inverse of the timescale of 
star formation, namely the timescale necessary to con- 
sume all the gas in the star forming region. More pre- 
cisely, T can be expressed as vri{t) where v is the real 
star formation efficiency, and r\(t) is a parameter taking 
into account the stochastic nature of the star formation 
processes (Gerola et al. 1980) as shown in Matteucci and 
Chiosi (1983), where a detailed description can be found. 

The rate of gas loss via galactic winds for each element 
is assumed to be simply proportional to the amount of gas 
present at the time t: 



G iw = Wi G(t) X iw (t) 



(3) 



where X iw (t) = Xi{t) is the abundance of the element i 
in the wind and in the interstellar medium (ISM); Wi is 
a free parameter describing the efficiency of the galactic 
wind and expressed in Gyr~ x . We have considered normal 
and differential winds: in the first case, all the elements are 
lost at the same efficiency (wi = w). In the second case, 
the value of Wi has been assumed to be different for dif- 
ferent elements: in particular, the assumption has been 
made that only the elements produced by type II SNe 
(mostly a-elements and one third of the iron) can escape 
the star forming region. We have made this choice follow- 
ing the conclusions of Marconi et al. (1994), who showed 
that models with differential winds can better explain the 
observational constraints of blue compact galaxies in gen- 
eral. The conditions for the onset of the galactic wind are 
studied in detail and described in paragraph 3.2. 

The chemical evolution equations include also an ac- 
cretion term: 



(Gi)inf = C 



(l.) m /e- (t/T) 

Mr 



where (Xi)i n f is the abundance of the element i in the 
infalling gas, assumed to be primordial, r is the time scale 
of mass accretion to is the galactic lifetime, and C is a 
constant obtained by imposing to reproduce Ml at the 
present time to- The parameter r has been assumed to be 
the same for all dwarf irregulars and short enough to avoid 
unlikely high infall rates at the present time (t = 0.5 • 10 9 
years) . 

Finally, the initial mass function (IMF) by mass, <j)(m), 
is expressed as a power law: 



4>{m) = Am- {1+X ^ 



(5) 



We considered two cases for the IMF: 

1. the exponent x — 1.35 over the mass range (0.1 -j- 
100) M Q (Salpeter 1955 IMF), and 

2. x — 1.35 over the mass range (0.1 -j- 2)M Q , and 
x = 1.70 over the mass range (2 4- 100)M o (Scalo 1986 
IMF). 

A is the normalisation constant, obtained with the follow- 
ing condition: 



(■mu 

A I mT x dm = l 



J m 



(6) 



3.1. Nucleosynthesis Prescriptions 



(4) 



The Ri{t) term of the equation (1) represents the stellar 
contribution to the enrichment of the ISM, i.e. the rate at 
which the element i is restored into the ISM from a stellar 
generation (see Marconi et al. 1994). This term contain all 
the nucleosynthesis prescriptions: 

1. for low and intermediate mass stars (0.8M Q < M < 
M up ) we have used Renzini and Voli's (1981) nucle- 
osynthesis calculations for a value of the mass loss 
parameter rj = 0.33 (Reimers 1975), and the mixing 
length any = 1-5. The standard value for M up is 8 
M ; 

2. for massive stars (M > 8M ) we have used Woosley's 
(1987) nucleosynthesis computations but adopting the 
relationship between the initial mass M and the He- 
core mass Mfle, from Maeder (1989). It is worth noting 
that the adopted M(Mn e ) relationship does not sub- 
stantially differ from the original relationship given by 
Arnett (1978) and from the new one by Maeder (1992) 
based on models with overshooting and Z = 0.001. 
These new models show instead a very different be- 
haviour of M(Mue) for stars more massive than 25 
Mq and Z = 0.02, but our galaxies never reach such 
a high metallicity; 

3. for the explosive nucleosynthesis products, we have 
adopted the prescriptions by Nomoto et al. (1986) and 
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Thielmann ct al. (1993), model W7, for type la SNe, 
which we assume to originate from C-0 white dwarfs 
in binary systems (see again Marconi et al. 1994 for 
details). 

As already said, nitrogen is a key element to understand 
the evolution of galaxies with few star forming events since 
it needs relatively long timescales as well as a relatively 
high underlying metallicity to be produced. The reason is 
that N is believed to be mostly a secondary element (sec- 
ondary elements are those synthesized from metals origi- 
nally present in the star and not produced in situ, while 
primary elements are those synthesized directly from H 
and He). In fact, N is secondary in massive stars, and 
mostly secondary and probably partly primary in low and 
intermediate mass stars (Renzini and Voli 1981). How- 
ever, some doubts exist at the moment on the amount of 
primary nitrogen which can be produced in intermediate 
mass stars due to the uncertanties related to the occur- 
rence of the third dredge-up in asymptotic giant branch 
stars (AGB). In fact, if Blocker and Schoenberner (1991) 
calculations are correct, the third dredge-up in massive 
AGB stars should not occur and therefore the amount of 
primary N produced in AGB stars should be strongly re- 
duced (Renzini, private communication). 

As a consequence, the only way left to produce a rea- 
sonable quantity of N during a short burst (no longer 
than 20 Myr) is to require that massive stars produce a 
substantial amount of primary nitrogen. This claim was 
already made by Matteucci (1987) in order to explain 
the [N/O] abundances in the solar neighbourhood. Re- 
cently, Woosley et al. (1994) have shown that massive stars 
can indeed produce primary nitrogen, and Marconi et al. 
(1994) and Kunth et al. (1995) have taken into account 
this possibility. In this paper we also have considered this 
possibility and, since the Woosley et al. data are not yet 
available we have used the same parametrization adopted 
by Marconi et al. (1994) and Kunth et al. (1995). 



3.2. Energetics 

Our model presents a new formulation of galactic winds 
relative to previous published models of this type (Mat- 
teucci and Tosi 1985 and Marconi et al. 1994). In partic- 
ular, we adopted the prescriptions developed for ellipti- 
cal galaxies (Matteucci and Tornambe 1987; Gibson 1994, 
1997): galaxies develop galactic winds when the gas ther- 
mal energy Eg h (t), exceeds its binding energy E b g (t), i.e. 
when: 



A detailed treatment of the energetics of the ISM is 
considered, in order to compute the gas thermal energy: 

Ef{t) = Ef N {t)+E? w {t) (8) 

we calculated the energy fraction deposited in the gas by 
stellar winds from massive stars El^t), and by super- 
nova explosions Eghty). Here the supernova contribution 
is given by: 



where: 

E^nii — / e SNRsNii(t') dt' 
Jo 



^SNIa 



Ef w {t) 



f 

Jo 



esNRsNia(t') dt' 



<f)(m)tp(t')e w dm dt' 



(9) 



(10) 



(11) 



(12) 



Jl2 



RsNii(t) and RsNia(t) are the rates of supernova (II and 
la) explosion, (j){m) is the IMF and ip(t) is the star for- 
mation rate. The type la and II SN rates are calculated 
according to Matteucci and Greggio (1986). The energy 
injected and effectively thermalized into the ISM from su- 
pernova explosions (esjv) and stellar winds from massive 
stars (e w ), are given by the following formulas: 



£SN = VSN E 



rjw E u 



(13) 



(14) 



Ef(t) > E b Jt) 



(7) 



where E is the total energy released by a supernova ex- 
plosion, Eyj is the energy injected into the ISM by a typical 
massive star through stellar winds during all its lifetime, 
and rjs n and r] w are the efficiencies of energy transfer from 
supernova and stellar winds into the ISM, respectively. We 
adopted a formulation for €sn and e w which is described 
in detail in appendix. 

Our formulation refers to an ideal case characterized by 
an uniform ISM, and no interaction with other supernova 
explosions or interstellar clouds. When a supernova ex- 
plodes the stellar material is violently ejected into the ISM 
and its expansion is gradually slowed down by the ISM. 
The energy released by the supernova explosion is assumed 
to be Eq — 10 51 erg, while the energy effectively trans- 
ferred and thermalized into the ISM is given by eq. (13) 
and depends on the assumed value for t]sn- The main pa- 
rameters used to derive 77s at, as described in the appendix, 
are: the initial blast wave energy E 5 i — i? /(10 51 erg), the 
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interstellar gas density no, and the isothermal sound speed 
in the ISM c 0:6 = c /(10 6 cm s _1 ). 

For typical values of these parameters E51 = 1, no = 
1 cm -3 and cq.q ~ 1 we obtained the following range of 
possible values for t]sn' 



VSN 



0.13-^0.007 



(15) 



In particular, we have adopted an intermediate value of 
Vsn — 0.03 which we assume to be the typical SN energy 
transfer efficiency. 

The energy injected into the ISM by a typical mas- 
sive star through stellar winds during all her lifetime is 
estimated to be: 



E w — L w tMS 



(16) 



where tms is the time the star spends in Main Sequence 
and L w is the stellar wind luminosity. By using the argu- 
ments developed in the appendix, we obtained: 



no/r \2/3 — 1/2 —5/2 

r\ w = 0.3(L 36 ) Z/J n CQg' 



(17) 



if we assume a 20M Q star as a typical massive 
star contributing to the stellar winds then L^q = 
L w /(10 36 ergs- 1 ) = 0.03, t M s = 7.9 • 10 6 yr with n = 
1 cm~ 3 and c , e = 1, we obtain: 

E w ~ 10 49 erg and rj w = 0.03 (18) 

Therefore, in our formulation we adopted: 

e S N = 0.03 • 10 51 erg (19) 

and 



0.03 • 10 49 erg 



(20) 



It is worth noting that the assumed value for E w is in 
agreement with that calculated by Gibson (1994) for a 
star of initial mass M — 20M Q . 

The introduction of dark matter halos with variable 
amounts and concentrations of dark matter is considered 
when we compute the binding energy of interstellar gas, 



E b = W L (t) + W LD (t) 



(21) 



The two terms on the right of the equation take in account 
the gravitational interaction between the gas mass M g (t), 
and the total luminous mass of the galaxy Mj,(i), and 
between the gas mass and the dark matter M<j: 



W L (t) 



-0.5 G 



M g (t)M L (t) 



(22) 



W LD (t) = -G WLD 



M g {t)M d 



(23) 



where w LD ~ ±S[1 + 1.375]. 



G is the gravitational constant, and S — r^ /r<j is the ratio 
between the effective radius of luminous matter and the 
effective radius of dark matter. These equations are taken 
from Bertin et al. (1992) and are valid for S defined in the 
range 0.10 -j- 0.45. 

It is worth noting that the original formulation of Bertin et 
al. (1992) was thought for massive elliptical galaxies, and 
that it is not necessarily the right one for dwarf irregulars. 
However, we used such a formulation since theoretical for- 
mulations for the binding energy of dwarf irregulars are 
not available. 

4. Model results 

Our study of the evolution of BCG and DIG starts from 
the results obtained by Marconi et al. (1994): i) the star 
formation is assumed to proceed in short and intense 
bursts of activity, which may induce galactic winds; ii) 
metal enriched and normal winds have been considered. In 
particular, for the metal enriched winds the assumption is 
made that galactic winds carry mostly the nucleosynthesis 
products of supernovae of type II, with the consequence of 
removing elements such as oxygen but not elements such 
as helium and nitrogen which are mainly produced in low 
and intermediate mass stars, and ejected through stellar 
winds. 

The novelty is that i) galactic winds are powered both 
by supernova explosions (SNII and SNIa) and stellar winds 
from massive stars, and that ii) we consider the presence 
of dark matter halos in these galaxies. 

In order to understand the observed distribution of 
N/O, C/O versus O/H, and of [O/Fe] versus [Fe/H], we 
computed different galaxy models by varying some param- 
eters such as the number of bursts, the duration of each 
burst, the star formation efficiency, the galactic wind effi- 
ciency, the dark matter mass and distribution, and finally, 
the IMF exponent. 

We first considered three sets of standard models char- 
acterized by 1, 3, 5, 7 and 10 bursts of a duration of 20, 60 
and 100 Myr, whereas the other parameters were fixed by 
the results presented in Marconi et al. (1994). In partic- 
ular, the IMF is the Salpeter one, and the star formation 
efficiency is T — 1 Gyr^ 1 . The galactic wind is differential, 
namely Wi is different from zero only for the a-elements 
ejected through type II supernova explosions. In particu- 
lar, for the elements studied here, Wi is zero for N and 
He, whereas it is Wi — 10 for O, Nc, Na, Mg and Si, is 
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Wi = 0.97 • 10 for C, since carbon is partly produced and 
ejected by low and intermediate mass stars through stel- 
lar winds, and is Wi = 0.3 • 10 for Fe since only ~ 30% 
of this element is produced by SNe II. The ratio between 
dark and luminous matter is assumed to be 10 and that 
between the luminous and dark core radius is S — 0.3 for 
the three sets of standard models. Normal galactic winds 
instead require Wi = to be the same for all the chemical 
elements. From now on we will indicate the models with 
differential winds by a capital D, (e.g. Wi = \0DGyr~ 1 
indicating the value for the a-elements). 



Salpeter IMF, R=10, w,= l, S = 0.3, n b -2 



1=0,1 

r=o.5 
r=7 



Salpeter IMF, R=10, w,= l, 5 = 0,3, n b = 2 



r=o.i 

r=i 

r=7 



7 8 
12 + log(0/H) 



Fig. 2. Two 0.06 Gyr bursts models: Salpeter IMF, R = 10, 
S — 0.3, Wi — lDGyr^ 1 . We present the observed distribution 
of log(C /O) vs 12 + log(0/H) and the variation range of the 
star formation efficiency of our models: F = O.I-t-7 Gyr~ x . The 
best fit however is given by T = 0.1 -f- 1 Gyr^ 1 . 



6 7 8 9 

12 + log(0/H) 

Fig. 1. Two 0.06 Gyr bursts models: Salpeter IMF, R = 10, 
S = 0.3, Wi = lDGyr^ 1 . We present the observed distribution 
of log(N/0) vs 12 + log(0/H) and the variation range of the 
star formation efficiency of our models: T = 0.14 7 Gyr^ 1 . The 
best fit however is given by T — 0.5 -j- 7 Gyr -1 . 

Among our standard models, those with burst dura- 
tion Atb = 20 Myr never developed galactic winds, and 
those with Atb = 100 Myr predicted the destruction of 
the galaxy (i.e. when the thermal energy of gas is larger 
than the binding energy of the galaxy) when a large num- 
ber of bursts was assumed (in particular when n b > 7); 
on the contrary models with = 60 Myr never pre- 
dicted galaxy destruction and developed galactic winds 
for rib > 5. As a consequence of this, we chose to examine 
only models with Atb = 60 Myr. 

We started then from our standard models and varied 
several parameters such as the star formation efficiency, 
the wind efficiency, the amount and distribution of dark 
matter. We found that in order to reproduce the observed 
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-0,2 - Salpeter IMF, R=10, w,-l, S-0.3, n.,-2 

r=o,i 

r=o,3 

-0.4 - r=3 

r=7 



-3 -2 -1 

[Fe/H] 

Fig. 3. Two 0.06 Gyr bursts models: Salpeter IMF, R = 10, 
S — 0.3, Wi — \DGyr~ 1 . We present the observed distribution 
of [O/Fe] vs [Fe/H] and the variation range of the star forma- 
tion efficiency of our models: r = 0.1 -=- 7 Gyr" 1 . The best fit 
however is given by T = 0.3 4- 3Gt/r _1 . 
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spread of log(N/0) and log{C/0) versus 12 + log(0/H) 
and [O/Fe] versus [Fe/H], the most important parame- 
ter is the star formation efficiency. We computed two series 
of models characterized by two different IMF: the Salpeter 
(1955) and the Scalo (1986) function. To reproduce the 
observed spread in all three diagrams, we needed a star 
formation efficiency (r) between 0.1 and 7 Gyr -1 when 
using the Salpeter IMF, and between 1 and 10 Gyr^ 1 
when using the Scalo IMF. The best value for the wind 
parameter was Wi = \DGyr~ 1 . 

Consider first the Salpeter IMF case: we computed 
models with 1, 2, 3, 5 and 7 bursts, and tested the im- 
portance of the presence of the dark matter. Greater is 
the amount of dark matter, deeper is the galaxy potential 
and higher the gas binding energy. Therefore, the thermal 
energy required to develop the galactic wind has also to be 
higher. We found that if the amount of dark matter is the 
same as that of the luminous matter (R = Md/Mr, = 1) 
any model predicts galaxy destruction during the first 
starburst if the star formation efficiency is high (r > 5). 
This upper limit of T depends of course on the num- 
ber of bursts: models with many bursts, are character- 
ized by a great consume of gas through star formation 
activity, so they must have lower star formation rates to 
keep bound, than models with just one or two bursts. We 
previously said that we need T = 0.1 7 Gyr^ 1 to ex- 
plain the observed spread of the abundances ratios. So, 
let's consider the presence of a greater amount of dark 
matter: ten times the luminous matter (i? = 10). In this 
case none of the models characterized by 1 or 2 bursts and 
r = 0.1-j-7 Gyr~ x predicts galaxy destruction during star- 
bursts, whereas those with 3 or more bursts sometimes do. 
Considering an even greater amount of dark matter, fifty 
times the luminous matter (R = 50), none of the mod- 
els predicts galaxy destruction for any value of the star 
formation efficiency F = 0.1 -j- 7 Gyr^ 1 . 

In table [t] we summarize the range of variation of the 
star formation efficiency and the occurrence of the galac- 
tic wind for different values of the number of bursts and 
different amounts of dark matter for stable models (i.e. 
the galaxy never blows up). In column 1 is presented the 
number of bursts (n& = 1, 2, 3, 5, 7), in column 2 the 
range of variation of the star formation efficiency (r) and 
in column 3 we indicate if models develop the galactic 
wind or not. This same scheme is reproduced three times, 
for R = 1, 10 and 50. 

Therefore, one of the main results of our study is the 
fact that the presence of dark matter halos around dwarf 
irregulars and blue compact galaxies is required in order 
to avoid total destruction due to the energy injected by 
supernova explosions and stellar winds during starbursts. 



Table 1. Salpeter IMF models characterized by S = 0.3, 
u>i — \DGyr~ 1 : we indicate the range of variation of the star 
formation efficiency (F, expressed in Gyr~ ) for different val- 
ues of R and rib, and the development of the galactic wind 
(GW). 





R = 1 


R = 10 


R = 50 


rib 


r 


GW 


r 


GW 


r 


GW 


1 


0.1^5 


yes 


0.1^7 


no 


0.1 - 


-7 


no 


2 


0.1 ^3 


yes 


0.1^7 


yes 


0.1 - 


-7 


no 


3 


0.1 ^ 1 


yes 


0.1^5 


yes 


0.1 - 


-7 


no 


5 


0.1 -f- 1 


yes 


0.1 ^3 


yes 


0.1 - 


-7 


yes 


7 


< 0.1 


yes 


0.1 1 


yes 


0.1 - 


-7 


yes 



Moreover, the distribution of dark matter in our formu- 
lation, is described by the parameter S = r^/rd, and we 
find that S can vary between 0.1 and 0.4. 

Figure 1 shows the range of variation of the star forma- 
tion efficiency for standard models characterized by two 
0.06 Gyr bursts (occurring at t = 1, and 5 Gyr), R — 10, 
S = 0.3, Wi = lDGyr- 1 and the Salpeter IMF. We no- 
tice the peculiar sawtooth behaviour, indicating the al- 
ternation of the active star formation periods (during the 
bursts) and the quiescent periods (during the interbursts). 
During the starbursts, the star formation rate is high, and 
the SNe of type II dominate in the chemical enrichcment 
of O, Fe, N and C, while during the interbursts, when the 
star formation is not acting, only the elements like C, N, 
or Fe, are produced by low and intermediate mass stars. 
In this case the abundance of these elements increases rel- 
ative to the oxygen abundance, and this is exactly what 
we observe in all the figures presented here. The observed 
spread of N/O abundance ratio, reported in figure 1, is 
quite well reproduced if T = 0.5 -j- 7 Gyr -1 . 

Figures 2 and 3 show the C/O and the [O/Fe] abun- 
dance ratios relative to the same models, respectively. In 
these cases we can notice how the observed spreads are 
quite well reproduced if T = 0.1 -j- 1 Gyr^ 1 for C/O and 
if T = 0.3 -f- 3 Gyr- 1 for [O/Fc]. These differences are 
probably due to the few data available in literature for 
C and Fe abundances and to uncertainties present in the 
nucleosynthesis calculations. 

On the other hand, figures 4, 5 and 6 show the same 
models but with the Scalo IMF, reproducing the N/O, 
C/O and [O/Fe] abundance ratios. In this case the best 
results are obtained with r = 1 -j- 10 Gyr^ 1 . However, 
while the [Fe/O] observed spread is quite well reproduced 
(figure 6), the C/O and in particular the N/O abundance 
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ratio are worsely reproduced than by the Salpeter IMF 
models (figures 2 and 1). 

On the basis of our study we can conclude that mod- 
els with the Salpeter IMF are favoured relative to those 
adopting the Scalo one, since these latter do not explain 
all the spread present in the N/O vs O/H diagram. The 
models with the Salpeter IMF require a star formation ef- 
ficiency T = 0.1 -T- 7 Gyr^ 1 and a dark to luminous matter 
ratio R = 1 50. 

In our analysis of the best model we have also con- 
sidered the relation between metallicity (Z) and the gas 
mass fraction (p = M gas /M t ot with M to t = Ml + Ma). 
Matteucci and Chiosi (1983) discussed the problem of 
reproducing the observed spread existing in the Z — fi 
diagram and suggested that different wind rates, differ- 
ent infall rates and different IMF from galaxy to galaxy 
could equally well explain the spread. However, they did 
not consider the role played by the dark matter. On the 
other hand, Kumai and Tosa (1992) suggested that the 
observed spread could be explained with the presence of 
variable amounts of dark matter: they proposed that the 
dark matter fraction should vary from galaxy to galaxy 
as f D = M d /M tot = 0.40 -e- 0.95. Our results agree with 
Kumai and Tosa (1992) and, as one can sec in figure 7, 
where it is shown that the dark matter can vary between 
1 and 50 times the luminous matter (R = 1 50) in order 
to reproduce the spread. This means that the parameter 
fn should vary in the range fo = 0.50 0.98. 

The second main result of our study concerns the en- 
ergetics of interstellar gas: the ISM receives energy from 
both supernova explosions and stellar winds from massive 
stars. In figure 8 we can notice the interstellar gas ther- 
mal energy relative to its binding energy and the total 
galaxy binding energy. When the thermal energy equates 
the gas binding energy, the galaxy develops the galac- 
tic wind. When a supernova explodes it injects almost 
10 51 erg into the ISM, but just some percents (see para- 
graph 3.) of this energy are transformed into thermal en- 
ergy of the gas. The question about stellar winds from 
massive stars is still under debate. As already discussed, 
we considered here a typical massive star (M ~ 2OM ) 
and assumed that it may inject into the ISM something 
like 0.03 • 10 49 erg through stellar winds during all of its 
life, and our results in this case suggest that the total 
thermal energy due to stellar winds from massive stars, 
is negligible if compared to the component due to super- 
novae of type II and la. Supernovae of type la also do 
not contribute significantly, as shown in Figure 9 where 
the different contributions to the thermal energy due to 
the SNcII, the SNela and the stellar winds are presented. 




Fig. 4. Two 0.06 Gyr bursts models: Scalo IMF, R = 10, 
S — 0.3, Wi — lDGyr~ . We present the observed distribution 
of log(NfO) vs 12 + log(0/H) and the variation range of the 
star formation efficiency of our models: F = 1 -f- 10 Gyr^ 1 . 



Scalo IMF, R= 10 w.= l, S = 0.H, n„ = 2 



r=i 
r=ic 




12 + Iog(0/H) 



Fig. 5. Two 0.06 Gyr bursts models: Scalo IMF, R = 10, 
S — 0.3, Wi — lDGyr^ 1 . We present the observed distribution 
of log(C/0) vs 12 + log(0/H) and the variation range of the 
star formation efficiency of our models: F = 1 -f- 10 Gyr -1 . 
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Scalo IMF, R=10, w,= l, S=0.3, n„=2 
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Fig. 6. Two 0.06 Gyr bursts models: Scalo IMF, R = 10, 
S = 0.3, Wi = \DGyr~ 1 . We present the observed distribu- 
tion of [O/Fe] vs [Fe/H] and the variation range of the star 
formation efficiency of our models: F = 1 -j- 10 Gyr -1 . 
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Fig. 7. Two 0.06 Gyr bursts models: Salpeter IMF, 
T = 5 Gyr~ x S = 0.3, vn = lDGyr' 1 . We present the ob- 
served distribution of Z versus Log/j, (fx = M ga s/Mtot). The 
models are characterized by a different amount of dark mat- 
ter: R=l, 10, 50. 



This would mean that both stellar winds and SNela play 
a negligible role in the evolution of these systems. 




5 10 

time (Gyr) 



Fig. 8. Three 0.06 Gyr bursts model, characterized by: 
Salpeter IMF, r = 5, R = 10, S = 0.3, ua = \DGyr~ 1 . 
The bursts occur at t = 1,5 and 11 Gyr. We present the galaxy 
binding energy (BEq), the interstellar gas binding energy (Eg) 
and the interstellar gas thermal energy (Eth)- The time of oc- 
curence of the galactic wind correspond to the time at which 

E t h = Eg. 

In agreement with the Marconi et al. (1994) results, 
our results also favour differential galactic winds, but we 
find that the wind efficiency parameter W has to be lower 
(W < 1 Gyr^ 1 ) than in Marconi et al. (1994). 

Finally in figure 10 is reported the characteristic be- 
haviour of type II supernova rates: each peak corresponds 
to a burst of star formation. In figure 11 instead typi- 
cal type la supernova rates are presented and we can no- 
tice how type la supernovae explode also during the inter- 
bursts periods. Considering Salpeter models characterized 
by R = 50, S = 0.3, and Wi = IDGyr^ 1 , we find that 
the present value of SNela rate varies between 0.32 • 10~ 9 
and 0.67 • 10 -4 yr^ 1 , depending on the values of both the 
star formation efficiency and the number of bursts. In par- 
ticular, in table we indicate the range of variation of 
the present value of type la SNe rate (RsNia) for differ- 
ent values of the number of bursts (n& = 1, 2, 3, 5, 7) and 
r = 0.1 -j- 7 Gyr- 1 . 
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Fig. 9. Three 0.06 Gyr bursts model, characterized by: 
Salpeter IMF, r = 5, R = 10, S = 0.3, w; = lDGyr" 1 . The 
bursts occur at t = 1,5 and 11 Gyr. We present the SNela 
and the stellar winds contributions to the total thermal energy 
relative to the SNell contribution: note that the stellar winds 
contribution is smaller than the one due to the SNela. 



time (Gyr) 



Fig. 10. Three 0.06 Gyr bursts model, characterized by: 
Salpeter IMF, r = 5, R = 10, S = 0.3, w; = lDGyr' 1 . The 
bursts occur at t = 1,5 and 11 Gyr. We present the rate of 
type II supernova explosion as a function of time. The units of 
the SN rate are SNeyr^ 1 



Fig. 11. Three 0.06 Gyr bursts model, characterized by: 
Salpeter IMF, T = 5, R = 10, S = 0.3, Wi = lDGyr' 1 . The 
bursts occur at t = 1, 5 and 11 Gyr. We present the rate of 
type la supernova explosion as a function of time. The units of 
the SN rate are SNeyr' 1 



Table 2. Range of variation of type la supernova rates (Rs Nia) 
for different values of the number of bursts (rib). The star for- 
mation efficiency is T = 0.1 ~ 7 Gyr^ 1 . 



rib 


RsNia (yr 1 ) 


1 


0.32 • 10" 9 - 


-0.80 ■ 10" 5 


2 


0.37- 10" 6 - 


-0.31 ■ 10" 5 


3 


0.24 ■ 10" 5 - 


-0.82 ■ 10" 4 


5 


0.25 • 10" 5 - 


-0.68 ■ 10" 4 


7 


0.30 • 10" 5 - 


-0.67- 10" 4 



5. Summary 

In this paper we attempted to model the chemical evolu- 
tion of dwarf irregulars and blue compact galaxies. In par- 
ticular, we considered the presence of dark matter halos 
and of the different contributions to the interstellar gas 
energy due to supernovae (II and la) and stellar winds 
from massive stars. The comparison with the available 
data on the abundances and abundance ratios of elements 
such as He, N, C, O, and Fe allowed us to conclude that 
DIG and BCG must contain a substantial amount of dark 
matter in order to be gravitationally bound against the 
intense starbursts, that galactic winds powered by super- 
novae and stellar winds from massive stars are preferably 
enriched (differential), and that the energy contribution 
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from type la supernovae to the total thermal energy of 
the gas is negligible relative to the type II supernova com- 
ponent. We found also that the stellar wind contribution 
is negligible relative to that of the type II supernovae. 
In summary, our models suggest: 

— a number of bursts n& = 1 10; 

— a star formation efficiency T = 0.1 -j- 7 Gyr^ 1 ; 

— differential galactic winds mostly powered by super- 
novae of type II, with a wind efficiency parameter 
W < IGyr^ 1 ; 

— a ratio between dark matter and luminous matter R = 
l-j-50; 

— a dark matter distribution given by the ratio between 
luminous effective radius and dark effective radius S = 
0.1 0.4; 

— Salpeter IMF is favoured relative to Scalo IMF: in fact 
the observed abundance ratios are better reproduced 
by Salpeter models than by Scalo models, in particular 
the spread of N/O abundance ratio. 



the hot bubble has a very long cooling time and cools adi- 
abatically as it expands pushing the outer shell. We thus 
assume that its thermal energy varies as E t h oc R~ 2 and 
neglect a correction oc i~ 4 / 7 due to the effect of cooling 
(Cioffi, Mckee & Bcrtschinger, 1988). This is compatible 
with the other simplifying assumptions made in our evalu- 
ation of the SN efficiency. The remnant therefore expands 
as R s oc t 2 / 7 in this 'snow plow' phase. At a later time the 
pressure inside the cavity becomes equal to the external 
pressure and the remnant stalls. We define the stalling ra- 
dius as the radius at which the shell velocity is equal to 
co. As long as the SNR expands supersonically, the ISM 
outside the remnant is not aware of its presence and the 
fraction of the SN explosion energy which is not cooled off 
remains locked inside the cavity. At merger, the shell dissi- 
pates its remaining kinetic energy and its thermal energy 
(we assume that the shell temperature does not become 
lower than the temperature of the unperturbed ISM, so 
that its thermal energy is three times the kinetic energy 
at the stalling radius) which add to the thermal content 
of the ISM. The efficiency is thus given by the ratio of this 
shell energy to the explosion energy. With the assumptions 
above it is found: 
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Appendix 

In order to calculate the energy gained by the ISM 
after an isolated SN explosion we consider the expansion 
of the SNR through an uniform medium of number den- 
sity no and isothermal sound speed co = 10 6 co,6 cm s~ 4 . 
Initially the temperature of the remnant is very high and 
the energy radiated away is negligible as compared to the 
explosion energy Eq\ this phase is called 'adiabatic' and is 
described by the Sedov solution, R s oc t 2/>5 . As the rem- 
nant expands its temperature decreases and eventually the 
radiative losses from the swept-up dense shell become rele- 
vant at the cooling time given by (cf. Cioffi & Shull 1991): 



E 

t c = 1.49 x 10 4 — r# 



3/14 



4/7, 



yr 



n ' 'C 5/14 

where Eq — 10 51 i?5i erg and ( is the metal abundance 
relative to the solar abundance. 

For t > t c most of the mass and the kinetic energy of 
the SNR is contained by the thin, cold, massive expanding 
shell, while the interior cavity is filled with hot dilute gas 
that contains most of the thermal energy. At this stage 



0.12E 2 5 { 35 n- 4/35 C 8/35 4^ 



It may happen that the remnant pressure becomes lower 
than the external pressure before the merger. In this case 
the remnant enters in the 'momentum conserving' phase, 
Rs oc i 1 / 4 . Spitzer (1968) has shown that, in the limiting 
case in which the 'snow plow' phase is entirely skipped and 
the 'momentum conserving' phase starts just after t c , the 
efficiency is given by , where v c is the remnant velocity 
at t c . In our notations we get: 



VSN 



0.005K 



51 



C0,6 



Thus, the 'true' efficiency will be intermediate between the 
two above. Andersen & Burkert (submitted), in describing 
the evolution of the ISM of the dwarf galaxies, assume 
£sn = 0.08 independent of the environment, but find that 
the exact value is not of importance. 

Before to explode, stars earlier than B0 (M > 20 M ) 
suffer considerable mass loss (M <~ 10~ 5 — 10~ 6 M Q yr -1 ) 
through supersonic winds (v w ~ 10 5 — 10 6 km s _1 ), carv- 
ing large, hot bubbles (Weaver et al., 1977) and injecting 
large amounts of energy into the ISM. In fact, as the wind 
impinges on the ISM, a reflected shock forms which ther- 
malizes the wind itself. This rarefied, hot, pressurizzed gas 
fills the expanding cavity which is surrounded by a thin, 
cold, dense shell formed by the swept up gas of the ISM. As 
the expansion procedes the radiative losses increase and 



12 



the inner pressure decreases untill eventually the bubble 
stalls. As for the SNR case, we calculate the kinetic en- 
ergy of the dense shell pushed by the expanding bubble 
when it becomes sonic. To obtain the wind efficiency e w 
the shell energy is then divided by the total mechanical 
energy of the wind ^ms^w, where tMS is the time spent 
by the star on the main sequence, and L w = 10 36 L36 erg 
s _1 is the wind mechanical luminosity. From the classical 
'energy conserving' solution of Weaver et al. , R oc i 3 / 5 , we 
get: 

_ n o r 2/3 -1/2 -5/2 
Vw — U -'J- ,J 36 n c 0,6 ' 

To obtain the above relation we assumed Tms = 4-4L 36 
(McKee et al. 1984). As in the case of SNRs, the above 
efficiency represents an upper limit; in fact, 'momentum 
conserving' bubbles (R oc i 1 / 2 ) may be realized if a sub- 
stantial enhancement of radiative losses occurs for some 
reason such as cloudy medium (McKee et al. 1984) or non 
steady cooling (e.g. D'Ercole 1992). In this case the hot gas 
quickly cools and the shell is pushed by the ram pressure 
of the wind that impinges directly on it. The efficiency 
then becomes: 

r/ w — 4./ x iu ^ 36 n v 2Q0Q c 06 , 

where V2000 is the wind velocity in units of 2 x 10 8 cm s _1 
We finally point out that in principle the SN efficiency 
worked out above could not be valid for exploding stars 
more massive than 20 M Q because the remnant inter- 
acts with the bubble rather than with the uniform ISM. 
If however we consider 20 M Q as representative of mas- 
sive stars (L36 = 0.03), then the maximum bubble radius 
is rather small (~ 10 pc), shorter of the radius where 
the SNR becomes radiative (~ 15 pc) and much shorter 
than the radius at which it stalls (~ 60 pc). If, moreover, 
heat conduction is active, the bubble tends to be replen- 
ished during the red supergiant phase of the star (D'Ercole 
1992). We thus consider the above estimate of esN rather 
accurate for our purposes. 
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